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The ferromagnetism of the checkerboard-lattice Hubbard model at quarter filling is one of the few exact ferromagnetic 
ground states known in the family of Hubbard models. When the nearest neighbor hopping, t\, is negligible compared 
with the second neighbor one, ti, the system reduces to a collection of Hubbard chains. We find that the ID character is 
surprisingly robust as long as t t <h_. This phenomenon of dimensional reduction due to the geometrical frustration leads 
to peculiar magnetic orders with ID character for intermediate U and is responsible for odd-frequency superconducting 
states close to the magnetic boundary. 

KEYWORDS: checkerboard-lattice Hubbard model, flat-band, ferromagnetism, dimensional reduction, odd- 
frequency superconductivity 



Strongly correlated electron systems (SCES) with geomet- 
rical frustration show a variety of fascinating physics. 1_8) The 
Hubbard model on the two-dimensional pyrochlore lattice, 
called checkerboard-lattice (CB-lattice) Hubbard model, is 
one typical example and its intriguing properties, such as the 
metal-insulator transition at l/2-fiiling 5 ' 8) and the ferromag- 
netism at 1/4-filling 1,2 ' have been studied by many authors. 
The latter is known as the flat-band ferromagnetism. Mielke 
have rigorously proved that the complete ferromagnetic state 
is the ground state of the CB-lattice Hubbard model at 1 /4- 
filling at the symmetric point t\ = ?2, 1) where t\ and t<x are the 
nearest and next nearest neighbor hoppings, respectively [see 
Fig. 1 (a)]. According to the numerical calculation, this fer- 
romagnetic state is extended to the case that t\ > ?2- 2) The 
electronic states including the ferromagnetism for t\ < t%, 
however, have not been investigated so far. 

In the SCES, the electrons move trying to avoid each other 
in space and time due to the strong local Coulomb repul- 
sion and it is difficult for them to condense into the isotropic 
pairing state where the on-site and equal-time gap function 
A(r = 0, t — 0) is finite with r and t being the relative coordi- 
nate and time of the Cooper pairs. One possible way to stabi- 
lize superconductivity in SCES is to form spatially anisotropic 
Cooper pairs, i.e. A(r = 0, f) = and such pairing states are 
considered to be realized in most of the SCES. 9) An alterna- 
tive way is to make the equal-time gap function A(r, t — 0) 
vanish. In other words, the gap function has odd-frequency 
dependence A(r, to) = -A(r, -co). The odd-frequency pair- 
ing state was first proposed by Berezinskii in the context 
of the 3 He. 10) After his proposal, theoretical studies on the 
odd-frequency pairing have been performed by several au- 
thors 1 1-1 8) and several models have been proposed as play 
grounds for realization of the odd-frequency pairing. Based 
on these studies, it has been pointed out that there are several 
favorable conditions for the odd-frequency pairing: (1) strong 
retardation effects, (2) geometrical frustration and (3) one- 
dimensionality. (1) When the retardation effects are strong, 
the effective interaction for the odd-frequency pairing is at- 
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Fig. 1. (Color online) (a) Schematic of the CB-lattice. The solid and 
dashed lines represent t\ and ti hoppings, respectively. The shaded region 
is the unit cell including two distinct sites A and B denoted by the filled and 
open circles, (b) Fermi surfaces for ?2 = 1.0 and ?i/?2 = 0.5 (solid lines) and 
for t\ = 1.0 and tijt\ = 0.5 (dashed lines). 



tractive for a wide range in frequency space and can dom- 
inate over the one for the even-frequency pairing. 18 * Strong 
retardation is realized in a SCES near a quantum critical point 
(QCP) 15) and also in an electron-phonon system with anoma- 
lously soft phonons. 18) (2) It has been shown that a geomet- 
rical frustration tends to suppress even-frequency spin-singlet 
pairing correlation and as a result, enhance the odd-frequency 
pairing correlation. 14 ' 17) (3) In the quasi-one-dimensional sys- 
tem, it has been shown that the odd-frequency spin-singlet 
p-wave pairing is favored when the one-dimensionality is 
strong. 17 ' In this case, the nodes of the gap function for the 
p-wave pairing do not intersect the Fermi surfaces and the 
pairing interaction is effective in a wide region of the mo- 
mentum space. Although it is obvious that the conditions (1) 
and (2) are satisfied in the CB-lattice Hubbard model near the 
QCP, we will show that the condition (3) is also satisfied in 
the model. 

In this letter, we investigate magnetic properties and super- 
conductivity of the CB-lattice Hubbard model at 1/4-filling 
by using the mean field approximation, exact diagonalization 

(FD^ method and random nhase annroxirnarion (RPA^ We 
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Fig. 2. (Color online) (a)-(c) qr-dependence of the largest eigenvalue of the 
bare susceptibility ^^(q) for t\ < t 2 at T = 0.02f 2 and (d)-(f) that for t\ > t 2 
at T = 0.02«i. 



show that the charge and magnetic ordered states with one- 
dimensional character are observed in the CB-lattice Hubbard 
model for t\ < t 2 . The spin fluctuations also have strong ID 
structure and promote the odd-frequency spin-singlet p-wave 
superconductivity. 

This remarkable property of the quasi-one-dimensionality 
is a consequence of the geometrical frustration of the present 
model, which may be characterized by dimensional reduc- 
tion. We will show that the odd-frequency spin-singlet p-wave 
state is stabilized in an extremely wide parameter range in this 
model. 

Let us first examine single-electron properties of the CB- 
lattice Hubbard model at 1/4 filling. When ti/t 2 = 0, there are 
two ID bands along k x and k y directions, each to be labeled 
by A- and B-band. Fermi surfaces for the A- and B-bands 
are respectively located along k x = ±3n/4 and k y = ±3n/4, 
since both t\ and t 2 are positive definite in this paper. Look- 
ing at the largest eigenvalue of the non-interacting paramag- 
netic susceptibility (denoted by XuJiQ) hereafter), the perfect 
nesting vectors of q = (n/2,q y ) and (q x ,n/2) give ridge-like 
structures, which are hallmarks of the one dimensionality. 
With increasing fi from zero toward t 2 , the A- and B-bands 
start to hybridize and we expect that the ID band structures 
would be modified significantly. Strikingly enough, though, 
the shape of Fermi surface does not change much, except for 
the small regions centered around the Fermi-surface crossing 
points at (k x ,k y ) = (±3^/4, ±3n/4) [see Fig. 1(b)]. As a re- 
sult, a major Fermi surface of a large hole pocket around T 
point looks square-like and retains one dimensionality even 
for t\jt 2 = 0.8. This is because the ID-like Fermi surface 
at 1/4 filling lies in the vicinity of the zone boundary of the 
first Brillouin zone where A-B band mixings vanish exactly 
from the symmetry. 19 * The ID nature of Fermi surface is also 
seen in xf^iQ) as ridge-like structures along q x - n/2 and 
q y = 7t/2 lines as shown in Figs. 2-(a) and (b). Along the 
ridge line, there is a weak structure and the most divergent 
mode is located at q = (n/2, n/2) for < ti/t 2 < 0.58 and 
at about q = (n/2,0) and (0,n/2) for 0.66 < h/t 2 < 0.91. 
In many cases, the geometrical frustration suppresses the q- 
dependence of the spin fluctuations, leading to suppression 
of the even-frequency anisotropic pairing. 14 ' 17 ' On the other 
hand, in the CB-lattice, it produces the above-mentioned ro- 



bust ID features of the single-electron properties and pro- 
motes the odd-frequency pairing as explained later. 

On the other hand for < t 2 /t \ < 0.9, the Fermi surface is a 
single rounded-square electron pocket [Fig. 1(b)] and^' t (q) 
has an incommensurate peak located at q = (n/2 + 6, n/2 + 6) 
with 6 ~ 0.1. Note that the ridge-like structures around 
q x = n/2 and q y = n/2 develop with increasing t 2 /t\, see Figs. 
2(d) and (e). As we approach the Mielke point by changing 
?2/?i or ti/t 2 to unity, due to the geometrically frustrated hop- 
ping processes, the lower band is flattened andxf^(q) become 
nearly-structureless with a weak peak located at q = (0, 0) for 
fi/f2 ~ 1 as shown in Figs. 2 (c) and (f). Near the Mielke 
point, therefore, strong electron correlation is intrinsic even 
when U is rather small compared with t\ and t 2 . 

After switching on the Coulomb interaction U, within the 
mean-field approximation, the peak modes in^j l(q) become 
leading magnetic instabilities from the paramagnetic phase 
and second-order phase transition takes place at a certain 
value of U c for a fixed t\ jt 2 . When t\ jt 2 = 1 , we know that the 
Mielke's ferromagnetism is realized irrespective of U, mean- 
ing that U c (t\/t 2 — 1) = 0. Exact diagonalization (ED) study 
for 16-site CB-lattice Hubbard model under anti-periodic 
boundary conditions shows that a critical U c from the para- 
magnetic (S = 0) to the ferromagnetic (S = 4) states grows 
rapidly by changing t\/t 2 or t 2 /t\ from unity, like U c (f 2 /t\ = 
0.8) = 4.0537fi or U c (h/t 2 = 0.8) = 4.5105f 2 . 2 ' 19) Since 
the flat-band ferromagnetism would be suppressed when the 
bands become dispersive, we consider the magnetic phase di- 
agram in U-t\ or U-tz plane by applying the mean-field ap- 
proximation. For this purpose, a unit cell is extended to the 

2 V2 x 2 a/2 ones shown in the inset of Fig. 3 and set up 
charge-modulated antiferromagnets (AFs) of line-, plaquette- 
, and checkerboard-type as well as usual paramagnet (PM) 
and ferromagnet (FM), for possible mean -field ground states. 
The line- and plaquette-AF have q - (n/2, n/2) spin struc- 
ture and are made from different arrangements of classical 
charge and spin ordered chains along the f 2 -bond directions 
at l/4-filling. In the case of tz/h < 1 and relatively large U, 
the checkerboard-type AF ground state is expected, since the 
charge modulation in this state are so distributed as to gain 
kinetic energy through t\ hoppings. These states are consis- 
tent with the short-range correlations calculated by ED. Fig. 

3 shows the magnetic mean-field phase diagram in U-t\ (left 
panel) or U-tz plane (right) as well as the RPA instabilities. 
The self-consistent energy of the different ordered states are 
typically calculated for a 32x32 number of enlarged unit cells 
at T = O.Olfi for t 2 /h < 1 or at T = 0Mt 2 for tjt 2 < 1. 19) 
From Fig. 3, one can see that Mielke's ferromagnetic phase 
is extended for large- U region and various AF phases emerge 
between the paramagnetic and ferromagnetic phases. When 
< h/t 2 < 0.2, the self-consistent energy of the line- and 
the plaquette-states are almost degenerate owing to the one- 
dimensionality of the system, which may lead to the discrep- 
ancy between the phase boundary estimated by the mean-field 
calculation at T = 0.01f 2 and that by RPA at T = 0.02f 2 . It 
should be also mentioned that, because of the size limitation 
of the unit cell in the mean-field approximation, we could not 
reproduce the second-order phase boundary from PM to the 
incommensurate state found in^^j for t\ > t 2 . It is natural 
to expect that there would be successive phase transitions be- 
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Fig. 3. (Color online) Mean-field phase diagram for (j < t\ (left-half panel) 
and t2 > t\ (right). The open circles and diamonds are, respectively, contin- 
uous and discontinuous phase-transition points within the present numerical 
accuracy. RPA-instability lines are also depicted by the dotted lines. Inset fig- 
ures (a)-(c) show the unit cells of various charge-ordered AF state, where the 
circular radius stands for the magnitude of charge density and each color zero 
(gray), positive (black), or negative (white) spin polarization. 



shown in Fig. 3. 

Now, let us turn our attention to the superconductivity. We 
investigate the superconductivity by solving the following lin- 
earized Eliashberg equation within the RPA, 



' k' a'P' 



k') 



X 



GZ(-k')A a ^(k')G^(k'), 



<0) , 



(1) 



where G^(k) is the noninteracting Green's function, &. a p(k) 
is the gap function , V a p(q) is the effective pairing interaction 
and A is the eigenvalue of the eigenvalue equation (1) which 
represents the strength of the superconducting correlation and 
reaches unity at T = T c . In the RPA, V a p(q) is given as, 



V{q)^r 1 U 2 x s {q)- l -U 2 x c {q) + U, 



(2) 



where 77 = 3/2(-l/2) for the spin-singlet (-triplet) state and 
the spin (charge) susceptibility x s(c) (l) is given as, %"*- c \q) = 
[l-(+)Ux (0) (q)V l X (0) (q) with the bare susceptibility x%(q) = 
~T/NZ k G^(k)G ( ®(k - q). We use the abbreviations k = 
(k,is n ) and q = (q,i(o m ), where s„ = (2n + 1)ttT and a>„, = 
2mnT. In the numerical calculations, we use the 128 x 128 
fc-meshes in the 1st Brillouin zone and 512 Matsubara fre- 
quencies (-51l7rr < e n < 51l7rr) and set T = 0.02fi and 
T - 0.02?2 for 1 1 > ?2 and t\ < ?2, respectively. 

Fig 4 shows the k- and the e„-dependence of the obtained 
gap functions for several typical parameter sets in the cases 
that t\ < ?2, ?i > ?2 and t\ ~ ?2, where the values of U are 
chosen as the largest eigenvalue A ~ 1. We note that in the 
present case, the effective pairing interaction given in eq. (2) 
can be approximated by V(q) ~ r]U 2 x s (q), that is, the super- 
conductivity is driven by the spin fluctuations because the spin 
fluctuations always dominate over the charge fluctuations and 
are strongly enhanced by the Coulomb interaction U when 
A ~ 1. With these in mind, we discuss the pairing symme- 
try and mechanism. It is found that the pairing symmetry for 
t\ = 1.0 and t2/t\ = 0.5 is even-frequency spin-singlet ex- 
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Fig. 4. (Color online) (a), (b) and (c) fc-dependence of the lower band di- 
agonal component of the gap function A(k,i£„) with the lowest Matsubara 
frequency e„ = nT for t\ = 1.0 and tilt\ = 0.5, t2 = 1.0 and f[/?2 = 0.5 and 
t2 = 1 .0 and t[/t2 = 0.9, respectively, The black and green (gray) lines denote 
the Fermi surfaces and the nodes of the gap function, respectively and the ar- 
rows schematically represent the typical pair scatterings, (d) £„ -dependence 
of that on the Fermi surface k = kf ~ (3n/4, 0) 



effective pairing interaction for the spin-singlet state given in 
eq. (2) has a sharp peak at q ~ (n/2, n/2) reflecting the struc- 
ture of^ (0, (^) shown in Fig. 2 (d), strong repulsive pair scat- 
terings with momentum transfer q ~ (n/2, n/2) take place 
[see Fig. 4 (a)]. Then, the gap function has different signs at 
the segments on the Fermi surface which are connected to 
each other by q ~ (n/2, n/2). 

On the other hand, the pairing symmetry for ?2 = 10 and 
fi/?2 = 0.5 is odd-frequency spin-singlet /?-wave (OSp) as 
shown in Figs. 4 (b) and (d). Since the effective pairing inter- 
action for the spin-singlet state given in eq. (2) has ridge-like 
structures along q x = n/2 and q y = n/2, strong repulsive pair 
scatterings with momentum transfer q ~ (n/2, q y ) and q ~ 
(q x , n/2) are important. Therefore, electrons on the Fermi sur- 
face with k x 37r/4 (k y 37r/4) are scattered to the other 

section of the Fermi surface with k x ~ 3n/4 (k y 37r/4) and 

vice versa [see Fig. 4 (b)]. Then, the gap function has sign 
change between the sections of the Fermi surface which are 
connected to each other by q ~ (n/2, q y ) and the resulting fc- 
dependence of the gap function is p v -wave. We note that the 
p x - and /^,-wave states are degenerate in the CB-lattice. Since 
the relation A(fc, is„) - A(-k, -is„) has to be hold for the 
spin-singlet state, the frequency-dependence of the gap func- 
tion is odd. 

In contrast to the above-mentioned two cases, near the 
Mielke point, the spin-triplet and spatially isotropic state, 
i.e., the odd-frequency spin-triplet s-wave (OTs) pairing is 
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point, the lower bandwidth is narrow and the q-dependence 
of x (?) is weak as shown in Figs. 2 (c) and (f). Then, the 
q-dependence of the effective pairing interaction for the spin- 
triplet state given in eq. (2) is also weak and V(q) < for any 
q which lead to strong attractive pair scatterings. Thus, the 
gap function shows no sign change and is of s-wave. Since 
A(fc, ie n ) = -A(-fc, -is„) has to be hold for the spin-triplet 
state, the frequency-dependence of the gap function is odd. 
Two distinct odd-frequency pairing states, the OSp and OTs 
states are realized due to the spin fluctuations with quite dif- 
ferent features characterized by the strong one-dimensionality 
on one hand and by the weak q-dependence on the other hand. 

Finally, we show the superconducting phase diagram in the 
t\- and t2-U planes in Fig. 5, where the superconducting phase 
boundary are defined as the points at which the largest eigen- 
value of eq. (1) A reaches unity. For t\ > *2, the ESs pairing is 
realized in the wide range of ?2/?i near the magnetic ordered 
phase with the ordering vector q ~ (jt/2,n/2) [see Figs. 2 

(d) and (e) and Fig. 3]. On the other hand, for t\ < ti, the 
OSp state is realized in the extremely wide range of ?i/?2- 
This situation is in a striking contrast to the usual quasi-one- 
dimensional Hubbard model, where the one-dimensionality 
is lost rapidly with increasing inter-chain hoppings leading to 
suppression of the OSp pairing correlation. The OSp pairing 
is realized also for t\ > t2 (0.73 < tz/ti < 0.83) because 
there the ID structure of£ (0) (g) develops as shown in Fig. 2 

(e) . Remarkably, the OSp pairing emerges already for mod- 
erate enhancement of the spin fluctuations, whereas the ESs 
pairing necessitates the very strong enhancement, which in- 
dicates the robustness of the OSp pairing. 17) It is worthwhile 
to note that for t\ < t%, the ESs state is second dominant in 
the wide range of fi/?2. The gap function for this ESs pair- 
ing, however, have no nodes on the Fermi surfaces in con- 
trast to the case that t\ > h because the Fermi surfaces are 
disconnected. 20 ' It is similar to the s ± -wave state discussed 
for the iron-based superconductors. 2 1_23) At 1/4-filling, the 
Fermi surface around k = (n, n) is smaller than that around 
k = (0,0) and the resulting T c for the ESs state is rela- 
tively low, while near 1 /2-filling, the Fermi surfaces around 
k = (n, ri) is comparable in size as that around k = (0, 0), 
and the resulting T c is relatively high (not shown). 19) Near the 
Mielke point, the OTs state is realized due to the spin fluc- 
tuations with featureless q-dependence. It should be pointed 
out that when the q-dependence of the spin fluctuations are 
weak, the mode-mode coupling effects of the various fluctu- 
ations are considered to become important. Thus, the vertex 
corrections as well as the self-energy corrections 3 ' neglected 
in the RPA may play significant roles for the stability of the 
OTs pairing near the Mielke point. We note that with decreas- 
ing temperature, the OTs region shrinks and the ESs region 
for t\ > t% and the OSp region for t\ < ti are extended toward 
the Mielke point because the q-dependence of the spin fluctu- 
ations gets stronger. Therefore, further investigations beyond 
the RPA study are needed to clarify whether the OTs pairing 
near the Mielke point eventually survives at lower tempera- 
tures. 

Here, we briefly comment on the re-entrant feature of the 
odd-frequency pairing previously reported. 15 ' 18) In the present 
study, we have not observed any clear re-entrant feature down 
to T = O.Olfi or T = O.Olfe- 
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Fig. 5. (Color online) Phase diagram on the t\- and ti-U planes, where the 
open circles, diamonds and triangles denote the ESs, OSp and OTs pairing 
instabilities, respectively. 



and the superconductivity in the CB-lattice Hubbard model 
at 1/4-filling with use of the mean field approximation, ED 
method and RPA. We have shown that the model exhibits the 
one-dimensional charge and magnetic orders, such as the pla- 
quette and line orders. The spin fluctuations have also the one- 
dimensional feature, i.e., the ridge-like structures in the mo- 
mentum space and drives the OSp superconductivity in the 
extremely wide range of t\lt%. These phenomena are due to 
the ID nature of the CB-lattice Hubbard model which is quite 
robust against the inter-chain hopping t\ in contrast to the case 
of the usual quasi-lD Hubbard model. 
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